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We study the Sivers effect in the transverse single spin asymmetries (SSA) for pion and kaon pro- 
I duction in semi-inclusive deep inelastic scattering (SIDIS) processes. We perform a fit of A^^''''' 

which, by including recent high statistics experimental data for pion and kaon production from HER- 
MES and COMPASS Collaborations, allows a new determination of the Sivers distribution functions 
for quarks and antiquarks with u, d and s flavours. Estimates for forthcoming SIDIS experiments 
at COMPASS and JLab are given. 
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I ■ I. INTRODUCTION 

&' 

\ In Refs. E, S], we studied the transverse single spin asymmetry A'^'^'^ observed by the HERMES Q and 
COMPASS [j| Collaborations in polarized SIDIS processes, ip{S) £' hX. The quality and amount of the data 
I allowed to perform a rather well constrained extraction of the Sivers distribution function 0] for u and d quarks, 
^ . assuming the existence of a symmetric and negligibly small Sivers sea. Similar analyses and extractions were performed 
by other groups 0, B ■ Although all these results were relevant and significant as the first determination of the 
, Sivers u and d functions, they were affected by the low statistics of the experimental data available at that time: in 

■ fact, COMPASS asymmetries were limited to charged hadron production, as no hadron separation was performed, 
CN ; while the HERMES data were only given for charged pion production. Recently, much higher statistics data on 
10 ' the A^j^'^'^'^ "^^-^ azimuthal asymmetries for SIDIS have become available: in Ref. [l^ the HERMES Collaboration 

, presents neutral pion and charged kaon azimuthal asymmetries, in addition to higher precision data on charged pion 

■ asymmetries; moreover, Refs. [ill [13 show the COMPASS Collaboration measurements for separated charged pion 
and kaon asymmetries, together with some data for Kg production. 

It is then timely and natural to reconsider the analysis performed in Ref. [2 in order to increase our understanding 
of the properties of the Sivers function. In particular, reduced error bars and hadron separation in both the HERMES 
and COMPASS sets of experimental data allow a better determination of the u and d flavour Sivers distribution 
' functions and, most importantly, a flrst insight into the sea and strange contributions to the Sivers functions, namely 
^^/iz/pT , A^/j/pT , A^/s/pT and A^/j/pt . 

Our strategy is the following. First we evaluate the impact of the new data with respect to the old data sets: as 
we shall explain in Section lllli using the same unpolarized fragmentation functions as in Refs. [ij, [2| would give a 
high quality fit as far as pion asymmetries are concerned, but fails to describe the kaon data. Instead, the use of a 
different, more recent set of fragmentation functions [l^, based on a global analysis of pion and kaon production, 
will prove to be crucial to reach a successful description of pion and kaon data simultaneously. With a simple ansatz 
parameterization of the Sivers functions, we will then perform a simultaneous fit of both HERMES and COMPASS 
data sets on A''^''^^^'^-^'^ for pion {tt^, 7r°) and production. Wc do not include in our fit the COMPASS data on 
Kg [13] as the corresponding fragmentation functions are not so well established and can be obtained from those for 

K^ only adopting further assumptions; wc shall rather estimate A^^'^'^'^ for Kg, and compare it with COMPASS 
data, using the Sivers functions obtained by fitting all other data sets, and assuming exact SU(2) invariance to derive 
the quark fragmentation functions into Kg. 

The above procedure will allow to determine the valence and sea proton Sivers functions, which will be used to 
provide estimates for the analogous single spin asymmetries that will soon be measured at JLab (operating on proton, 
neutron and deutcron targets) and at COMPASS (operating on a proton target). Notice that the JLab measurements 
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will provide vital information on the large x behaviour of the Sivcrs distribution functions, yet undetermined from 
present SIDIS experiments, as explained in Section HVl 



II. FORMALISM AND PARAMETERIZATION 

The SIDIS transverse single spin asymmetry (SSA) A'jl^^''^'^"'^^^ measured by HERMES and COMPASS is defined 
as (see Fig. [l]for the definition of the azimuthal angles) 

.sin(0fc-0s) _ 9 / d(j)s d(j)h [rfo-^ - da^] sinjcph - <j>s) x 
Jdcj,sdcj)h[daT +dai] ' ^' 

and shows the azimuthal modulation triggered by the correlation between the nucleon spin and the quark intrinsic 
transverse momentum. This effect is embodied in the Sivers distribution function A^fq/pT{x,k±), which appears in 
the number density of unpolarized quarks q with intrinsic transverse momentum k±^ inside a transversely polarized 
proton , with three- momentum P and spin polarization vector S, 

fg/pi{x,k±) = fg/p{x,k±_) + ^A^fg/pT{x,k±) S ■ {P X kj_) , (2) 

where fq/p{x, k±) is the unpolarized x and kj_ dependent parton distribution, and the mixed product S ■ {P x k±) 
explicitly gives the azimuthal dependence mentioned above. Notice that the Sivers function is also often denoted as 
f^j!{x, k±) [ij]; this notation is related to ours by [l^ 

A%/p, (x, k^) - fi^{x, k^) . (3) 

Tllp 

The "weighting" factor sm{(j)h — (j^s) in Eg . IT]) is appropriately chosen to single out, among the various azimuthal 
dependent terms appearing in [da^ ~ da^] |16l.ll7l|. only the contribution of the Sivers mechanism. By properlytaking 
into account all intrinsic motions this transverse single spin asymmetry can be written, at order {k±/Q), as [4] 

d(ps d(j)h d'^kj_ A^fq/pr{x,k±)sin{ip - (j>s) D'^{z,pA^)s\n{cj)h - (j)s) 

^UT — I. _T-_f„^fn ' 



d(j)s d(j)h d^k± fq/pix, k±) D'^{z,pa_) 



(j)s and (jjh are the azimuthal angles identifying the directions of the proton spin S and of the outgoing hadron 
h respectively, while if defines the direction of the incoming (and outgoing) quark transverse momentum, k± = 

fci(cos(/3,sin(/3, 0), as shown in Fig. [TJ '^""^^3 " is the unpolarized cross section for the elementary scattering iq iq, 

dQ^ s2 Q4 ' (5) 

where s, t = —Q^ and it are the partonic Mandclstam invariants. 

Finally, D^{z^p±) is the fragmentation function describing the hadronization of the final quark q into the detected 
hadron h with momentum Ph (see Fig. [Ij; h carries, with respect to the fragmenting quark, a light-cone momentum 
fraction z and a transverse momentum . 

In our analysis we shall consider u, d and s flavours for quarks and antiquarks. The Sivers function is parameterized 
in terms of the unpolarized distribution function, as in Ref. Q], in the following factorized form: 

A'^/g/pT (X, fci) = 2Afg{x) h{k^) fq/pix, k^) , (6) 

with 

Nq{x) = N, (1 - xf^ ^"^'^y. ' (7) 

aq^Pq 

Mfc^) = V2i|f e-'^l/^^?, (8) 

IVli 
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where Nq, aq, (3q and Mi (GeV/c) are free parameters to be determined by fitting the experimental data. Since 
h{k±) < 1 for any k± and |A/'q(a;)| < 1 for any x (notice that we allow the constant parameter Nq to vary only inside 
the range [—1,1]), the positivity bound for the Sivers function. 



2/q/p(a;,fc^) 



is automatically fulfilled. We adopt the usual (and convenient) Gaussian factorization for the unpolarized distribution 
and fragmentation functions: 



fq/p{x, k±) — fq{x) ^^^2 



and 



D'^{z,p^)^D'^{z)-^e-A/ipl) , (11) 

with the values of {kjj and {p\) fixed to the values found in Ref. [if by analysing the Cahn effect in unpolarized 
SIDIS: 

(fci) = 0.25 (GeV/c)2 {pD = 0.20 (GeV/c)^ . (12) 



Notice that the Gaussian distributions limit the effective action of intrinsic motion to fc^ < \/{k±) a nd p \ < ^/ (p^, ) , 
which is the region of validity of the TMD factorized expressions in Eq. (jl]), Pt ~ fc^ ~ Aqcd <C Q fla. [l9l. 1201 . 

The parton distribution functions (PDF) fq{x) and the fragmentation functions (FF) D^{z) also depend on via 
the usual QCD evolution, which will be taken into account, at leading order (LO), in all our computations. 

Before fitting the data on the Sivers asymmetries a few comments on the quark hadronization are necessary. While 
most of the available sets of fragmentation functions describe rather well the pion multiplicities observed at HERMES, 
many of them fail to reproduce the kaon multiplicities in SIDIS production. The main reason is the role of the strange 
quarks, which is often not well established: for example, one expects that mesons can be abundantly produced 
by s quarks, via creation from the vacuum of a light uu pair, rather than by u quarks, via creation from the vacuum 
of a heavier ss pair. Such a feature is particularly emphasized in the set recently obtained by de Florian, Sassot, 
Stratmann (DSS) [H, which has Df^{z) > D^^ {z) over the whole z range. This is shown in Fig. [2l where the LO 
DSS fragmentation functions (solid lines) are compared with those proposed by Kretzer (KRE) [21| (dashed lines) 
and by Hirai, Kumano, Nagai and Sudoh (HKNS) [22] (dotted lines). The DSS set, which is determined by fitting all 
presently available multiplicity measurements, both for pions and kaons, is indeed the most suitable for our purposes. 

This can also be seen in a more quantitative way. We know that Kretzer's and other commonly adopted sets 
of fragmentation functions are able to describe pion production data, as shown, for instance, in Fig. 4 of Ref. [31 ■ 
However, Fig. 13 of Ref. [l^ shows instead that Kretzer fragmentation functions fail to reproduce charged kaon SIDIS 
multiplicities, and might not be adequate to reconstruct transverse single spin asymmetries corresponding to kaon 
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FIG. 2; The LO unpolarized fragmentation functions for it, d and s flavours, at = 2.4 (GeV/c)^, as given by Kretzer [2l[ 
(dashed lines), by Hirai, Kumano, Nagai and Sudoh [S^l (dotted lines) and by de Florian, Sassot and Stratmann [l3| (solid 
lines) . We show the fragmentation functions for 7r+ in the upper panel, and for K'^ production in the lower panel. Notice that 
the fragmentation function describing the probability of a K'^ originating from an s quark is much larger than that associated 
to a K'^ from a u quark, over the whole z range. The two solid curves in the central upper plot correspond to zDJ — zD^ 
(upper line) and zD'^ — zD'^^ (lower line). 



production. In fact, by using the Kretzer set for our fit, we would not be able to describe the kaon asymmetry data: 
to be more precise, we would obtain /d.o.f. = x^o/ — ^ ^O'" pions but Xdo/ — ^ foi' kaon production asymmetries. 
Estimates for asymmetries were presented in Ref. Hand the inadequacy of the Kretzer fragmentation functions 
was pointed out in several talks (see, for example, Ref. [23l|). The same conclusion has been confirmed, very recently, 
in Ref. [H. 

Let us now turn to the experimental data on kaon and pion Sivers azimuthal asymmetries measured by the HERMES 
Collaboration [l0| . The single spin asymmetry corresponding to K'^ production is, as a matter of fact, much larger 
than the analogous asymmetry for 7r+. Although one could naively expect, on the basis of u quark dominance, that 
K'^ and tt"*" asymmetries should be roughly the same, the presence of a large Df FF can help to understand the 
"puzzle" of the asymmetry. Indeed, if a non-negligible s Sivers function exists, then its action combined with a 
large Df fragmentation function can give rise to a significant difference between K'^ and 7r+ Sivers asymmetries. 



III. FIT OF SIDIS DATA AND EXTRACTION OF SIVERS FUNCTIONS 



The recent SIDIS experimental data on Sivers asymmetries for pion and kaon production give us the opportunity 
to study sea-quark Sivers functions for u, d, s and s quarks. The u and d quark Sivers functions alone were already 
studied in Ref. ; in the present analysis we will be able to improve the extraction of these functions and to present 
first estimates of the sea-quark Sivers functions. 

In order to evaluate the significance of the sea-quark Sivers contributions we first perform a fit of the SIDIS data 
using flavour independent ratios of the sea-quark Sivers functions with the corresponding unpolarized PDFs: that is, 
for u, d, s and s flavours we attempt an "unbroken sea" ansatz: 

^^fq/p^ (a;, k±) = 2J\fsea{x) h{k±) fq/p{x, k±) , 
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-'sea 



where Nsea, ctsea, /9sea are the same for all sea quarks, q = u, d, s and s. 

As the SIDIS data from HERMES and COMPASS have a limited coverage in x, typically a; < 0.3 - 0.4, the 
experimental asymmetries we are fitting contain very little information on the large x tail of the Sivers functions. 
In fact, our previous analysis of Ref. [4I showed that the parameters /3„ and Pd as determined by MINUIT best fit 
procedure are affected by very large errors. Therefore, as a first attempt, we assume the same value of (3 for all Sivers 
functions, setting (i^ea = = Pd = fi- 

Thus for the "unbroken sea" ansatz we have 8 free parameters: 



* sea 



Nu Nd N,, 

(14) 

P Ml (GeV/c) . 

For the purposes of our fit, we use the unpolarized parton distribution functions fq{x,Q'^) as given in Ref. [25j 
(GRV98LO) and the fragmentation functions DI^{z,Q^) as given in Ref. [T3| (DSS) - all evolved to the appropriate 
values - with the additional k± Gaussian dependences of Eqs. (fTU|) - P^ . While the choice of the DSS fragmentation 
functions is the one which best describes the large asymmetries observed for if + , the usj3of different sets of distribution 
functions, including the most recent analysis of s quark distributions from HERMES ^261 , would not affect our results 
significantly. For the Sivers functions, we use the functional forms of Eqs. ([HI)-®. Notice that the (unknown) 
evolution of these functions is assumed to be the same as for the unpolarized PDFs, fq{x, Q^)- 

By fitting simultaneously pion and kaon production data from HERMES [3 and COMPASS ^ we obtain an 
acceptable overall description of the experimental data, with x^gj = 1.16. The new experimental data give clear 
indications of the need of a non negligible sea-quark Sivers function, with Ng^a ~ —0.13 ± 0.03 sensitively different 
from zero. Although the total Xdof definitely good, a more careful examination of the results shows that while we 
achieve a perfect description of the 7r+ production data at HERMES [l3|, with — 1 per data point, the description 
of kaon production data is rather poor, with ~ 3 per data point for production at HERMES [l3|. This 
indicates that the "unbroken sea" ansatz fails to reproduce the differences between pion and kaon production, and 
clearly suggests the need of a parameterization which should allow for a more structured flavour dependence of the 
sea-quark Sivers functions. 

Including four new functions in our analysis would result in a substantial growth of the number of parameters and 
would consequently limit the usefulness of our parameterization. To keep the number of parameters under control, 
we define a simple "broken sea" ansatz by introducing four free parameters, Nu, N^, Ns, and Ng which give different 
sizes to the sea-quark Sivers functions, while keeping the same functional forms (uu = = cts = ctg = asea and 
Psea = Pu = Pd = P) ■ For thc "brokcn sea" ansatz fit we then have 11 parameters: 



Nu 


Nd 


Ns 




Nd 


Ns 




ad 




P 


Ml 


(GeV/c) . 



(15) 

The results we obtain for these parameters by fitting simultaneously the four experimental data sets on A^'^'^'' 
corresponding to pion and kaon production at HERMES [l3l and COMPASS [ll[, are presented in Table [fl together 
with the corresponding errors, estimated according to the procedure outlined in Appendix [XI The fit performed under 
the "broken sea" ansatz shows a remarkable improvement, especially concerning the description of kaon data. We 
now obtain x^ = 1-20 per data point for K'^ production at HERMES [l3|, while for pions we have x^ = 0.94 per data 
point, and a total Xdof — 1-00. In Table II we show the x^ per data point for pion and kaon production at HERMES 
and COMPASS, both for the "unbroken sea" and "broken sea" ansatze and adopting thc Krctzer and DSS FF sets. 
Notice that these values refer to the asymmetries as a function of x. 

The quality of our results is shown in Figs. [3] and 3] where our best fit to the SSA is compared with the experimental 
data from Refs. [l3| and 11 1: the SSAs are plotted as a function of one variable at a time, either a; or z or Pt, while an 



integration over the other variables has been performed consistently with the cuts of the corresponding experiment. 

In order to check the dependence on the set of unpolarized PDF's adopted, we have also performed the fit by using 
the CTEQ6L and the MRSTOILO [111 sets; in both cases, the quality of the fit and the central-value results for 
the asymmetries are so similar to those obtained with the GRV98LO set that they would be hardly distinguishable 
in Figs. [3] and H 
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TABLE I: Best values of the free parameters for the "broken sea" ansatz, Eq. (|15fl . Notice that the statistical errors reported 
in this table are not the errors given by MINUIT. As the parameters are strongly correlated, we determine them according to 
the procedure explained in Appendix A. The significant fluctuations in our results are shown by the shaded areas in Figs. |3] 
and SI 



xVd.o.f. = 1.00 



Nu = 0.04t°:g Ns = -0.40l°:S iV, = ito.ooo, 

/? = 3.46l:^-^^ Ml = 0.34«-?« (GeV/c)2 



The shaded areas in Figs. [3] and |4] (and in all subsequent figures where they are shown) represent statistical 
uncertainties and correspond to a 95.45% Confidence Level (CL): they are determined according to the procedure 
described in Appendix [X] Notice that further uncertainties of theoretical nature, intrinsic to our phcnonienological 
approach, arc present and might widen the size of the statistical bands. However, these are very difficult to assess: it 
suffices to recall that our analysis is performed assuming a simple factorized fcj, dependence in Eqs. ([6]), (fTO|) and (fTTjl . 
that the actual evolution of the Sivers function is unknown and that uncertainties in the fragmentation functions 
have not been taken into account. The functional form for the x-dependence of the Sivers functions used in the fit, 
Eq. ([7]), is a simple one and more structured dependences might allow better fits, with the statistical shaded areas 
covering better the experimental errors of the data. At this stage, considering the available experimental information 
and the remaining theoretical issues to be clarified, we do not think that further refinements of our analysis would be 
meaningful. 

Notice that in Fig. 2] we also show the results for 7r° at COMPASS, for which no data is so far available, computed 
using our extracted Sivers functions as given in Table HI Similarly we have computed A^!^'-'^'' for K'g production 
at HERMES and COMPASS and show them respectively in Figs. [3] and g] As the K\ is an equal mixture of = ds 
and = ds, we have assumed isospin invariance, writing the Kg EEs in terms of the ones - which are taken 



TABLE II: P^r data point for pion and kaon production at HERMES and COMPASS, both for the "unbroken sea" and 
"broken sea" ansatze, and adopting the Kretzer and DSS FF sets, x dependent data, integrated over z and Pt, are considered 
here. 
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FIG. 3: The results obtained from our simultaneous fit of the SIDIS A' 



sin(</>)i-'^s) 
UT 



Sivers asymmetries (solid lines) are compared 
with HERMES experimental data [10| for pion and kaon production (left and right panel, respectively). The shaded area 
corresponds to the statistical uncertainty of the parameters, see Appendix|X]for further details. For completeness, we also show 
the K% asymmetry, not measured at HERMES, which is the result of a computation based on our extracted Sivers function 
and the assumed fragmentation functions of Eq. (|16p . 
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FIG. 4: The results obtained from our fit (solid lines) are compared with the COMPASS measurements of ^ 
pion (left panel) and kaon (right panel) production off a deuteron target. The shaded area corresponds to the statistical 
uncertainty of the parameters, as explained in Appendix]^ The 7r° asymmetry, not measured at COMPASS, is the result of 
a computation based on our extracted Sivers functions. Also the K% asymmetry, although compared with data is not a 
best fit, but the result of our computation, using the assumed fragmentation functions of Eq. (|16p . 
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Our computation of the Kg asymmetry at COMPASS can be compared with the available data [l2|, as shown in the 
upper-right plots of Fig. [H Notice that these curves, contrary to the others in the same figure, are not best fits, but 
a simple estimate, based on the extracted Sivers functions and the adopted fragmentation functions of Eq. p6p . 
In Fig. [51 our results, obtained using the kaon fragmentation functions as given by dc Florian et al. in Ref. (isj 
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FIG. 5: The results obtained from our fit using the kaon fragmentation functions as given by de Florian et al. in Ref. [l^l (sohd 
hnes) are compared with the results we would find by using the KRE [23|(dotted lines) and HKNS [l^] (dashed lines) sets of 
fragmentation functions. 



(solid lines), are compared with the best fit we would find by using the KRE [2l| (dotted lines) and HKNS [l^ (dashed 
lines) sets of fragmentation functions. It is clear that the use of the new - strange-quark sensitive - fragmentation 
functions yields a much better agreement with the experimental measurements of the SIDIS azimuthal asymmetries 
for kaon production. 

The Sivers functions generated by our best fit procedure arc presented, at the scale = 2.4 (GeV/c)^, in Fig. [6l 
where we plot, on the left panel, the first fe^ moment defined as 

A%1, ix) ^Jd'k^^^ A^/,/,r (X, k^) = -ft^'Hx) , (17) 

and, on the right panel, the fc^ dependence of A^/^/pT at a fixed value of a; = 0.1. The highest and lowest dashed 
lines show the positivity limits |A^/| = 2/. 

Our results both confirm previous conclusions on the u and d Sivers distributions and, despite the still large 
uncertainties (see Table |T|, offer some new clear information about the so far unknown sea-quark Sivers functions. 
Let us comment in detail: 

• The HERMES data on kaon asymmetries, surprisingly large for , cannot be explained without a sea-quark 
Sivers distribution. In particular, we definitely find 

A%pT > (18) 

and confirm the previous findings for valence flavours 0, 0, H, 1 

A^/„/pT > A^/,/pT < . (19) 

There are simple reasons for the above results. The Sivers distribution function for s quarks turns out to be 
definitely positive, due to the large positive value of A^^'^'^ '^^^ for + ; notice that the value of Ns saturates 
the positivity bound \Nq\ < 1. Similarly, the positive sign of A^/„/pT is, essentially, driven by the positive tt"^ 
and A'+ SSAs and the opposite sign of A^/^/pt by the small SSA measured by COMPASS on a deuteron target. 
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X (GeV) 

FIG. 6: The Sivers distribution functions for u, d and s flavours, at the scale = 2.4 (GeV/c)^, as determined by our 
simultaneous fit of HERMES and COMPASS data (see text for details). On the left panel, the first moment x A'^f^^\x), 
Eq. (|17p . is shown as a function of x for each flavour, as indicated. Similarly, on the right panel, the Sivers distribution 
X A'^f{x, k±) is shown as a function of k± at a fixed value of x for each flavour, as indicated. The highest and lowest dashed 
lines show the positivity limits |A^/| = 2/. 

The u and d Sivers functions are also predicted to be opposite in the large Nc limit [2^ and in chiral models 

M- 

• The Sivers functions for u, d and s quarks, instead, turn out to have much larger uncertainties; even the sign 
of the u and s Sivers functions is not fixed by available data, while A^f^/pt appears to be negative. This could 
be consistent with a positive contribution from u quarks, necessary to explain the large asymmetry, which 
is decreased, for , by a negative d contribution. One might expect correlated Sivers functions for s and s 
quarks: we have actually checked that choosing A^/^/p? — ±A''^/j/pT slightly worsens the Xdof (from 1 up to 
about 1.1), but still leads to a reasonable fit. 
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FIG. 7: The Sivers distribution functions for u and d flavours, at the scale = 2.4 (GeV/c)^, as determined by our present flt 
(sohd hnes) , are compared with those of our previous fit 0] of SIDIS data (dashed fines) , where tv'^ and kaon productions were 
not considered and only valence quark contributions were taken into account. This plot clearly shows that the Sivers functions 
previously found are consistent, within the statistical uncertainty bands, with the Sivers functions presently obtained. 



• We notice that the Burkardt sum rule [3l| 

a a 

where, from Eqs. and (fT7|) . 

(fcl) = 



dx 



dkj_k^A fa/pr{x,k^) 



(SxP) 



is almost saturated by u and d quarks alone at ~ 2.4 (GcV/c)^: 

{kl) + (fci> = -nni (MeV/c) (fc!) + (fci) + (fcl) + (fcl) = -14t^3 (MeV/c). 

The individual contributions for quarks are: 





- 96i^° 


(MeV/c) 


(ki) = 


-1131^ 


(MeV/c) 


(kl) 


9+24 
— ^-11 


(McV/c) 


(kl) = 


9Q+20 


(MeV/c) 


{kl) 


- _4+ii 


(MeV/c) 


{kl) = 


-1 7+30 
-^'-8 


(MeV/c) , 



(20) 



(21) 



(22) 



(23) 



thus leaving httle room for a gluon Sivers function, 

- 10 < (kl) < 48 (MeV/c) , (24) 

in agreement with other similar results [s^, . The statistical uncertainties in the values given above have 
been computed as explained at the end of Appendix A. 
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FIG. 8: Estimates for the single spin asymmetry ^^y''*'' for pion and kaon production off a hydrogen target, which will 
be measured by the COMPASS Collaboration. 



• In Fig.[7]we compare the u and d flavour Sivers distribution functions, at the scale = 2.4 (GeV/c)^, obtained 
in our present analysis with the u and d flavour Sivers functions wc had found from our previous fit [31 , where tt^ 
and kaon productions were not considered, the Kretzer fragmentation function set was used, and only valence 
quark contributions were taken into account in the polarized proton. This plot shows that the Sivers functions 
previously obtained are consistent, within the statistical uncertainty bands, with the Sivers functions presently 
found. 



IV. ESTIMATES FOR FORTHCOMING EXPERIMENTS 



Using the Sivers functions determined through our fit, wc can give estimates for other transverse single spin 
asymmetries A^'^'^'^~'^^'^ which will be measured in the near future. Fig. [5] shows the results we obtain for the 
COMPASS experiment operating with a hydrogen target, adopting the same experimental cuts which were used for 
the deuterium target (Eq. (71) of Ref. Q). 

Forthcoming measurements at the energies of 6 and 12 GcV are going to be performed at JLab, on transversely 
polarized proton, neutron and deuteron targets. The obtained data will be important for several reasons; they will 
cover a kincmatical region corresponding to large values of x, a region which is so far unexplored by other SIDIS 
measurements. In particular, a combined analysis of HERMES, COMPASS and JLab SIDIS data will allow a much 
better determination of the f3 parameters, which control the large x behaviour of the Sivers distribution functions. 
In addition, the combined analysis of proton and neutron target events will help flavour disentangling and a more 
precise determination of u and d quark contributions. 

Our estimates for the JLab SSAs, for pion and kaon production off proton, neutron and deuteron targets, at 12 
GeV, are presented in Figs. [5HTT1 At this energy relatively large values are available and the TMD factorization, 
valid for Pt ~ fcj^ ~ Aqcd <C Q, should hold. The adopted experimental cuts for a proton or a deuteron target are, 
in terms of the usual SIDIS variables, the following: 

0.3 < Zh< 0.8 
0.05 <x^ < 0.7 
1 < Q2 < 8 (GeV/c)2 
1-5 <Eh< 3.5 GeV , 

whereas for a neutron target they are: 

0.3 < Zh < 0.7 
0.34 < 2/ < 0.9 
W2 > 2.3 GeV^ , 

where, at order (kj_/Q), one has — x and — z; tl 



0.05 <Pt< 1.5 GeV/c 
0.25 < , < 0.85 
W^>A GeV^ 



0.05 < a;^ < 0.55 

Q2 > 1 (GeV/c)2 (26) 



he exact relationships can be found in Ref. [l|. 
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Notice that these estimates, while well constrained by the available SIDIS data at small x values, might be less 
stringent at large x values: for example, relaxing the assumption of a unique (3 value for all flavours, would only 
marginally affect our present fit of HERMES and COMPASS data, but would much widen the uncertainty band 
above x ~ 0.4. 

We have given complete sets of estimates for charged and neutral pions and kaons. While the tt^, 7r° and com- 
putations originate from a consistent procedure which uses fragmentation functions and Sivers distributions obtained 
by fitting data involving the same particles, the Kg estimates might be affected by a greater uncertainty about its 
fragmentation functions, which require the assumptions of Eq. (jl6p . This can be seen in the comparison between data 
and computations of Fig. [4l 

We have also computed estimates for JLab operating at 6 GeV, with the corresponding kinematical cuts: 

0.4 <zh< 0.7 
0.1 <x^ < 0.6 
Q2 > 1 (GeV/c)2 
l<Eh<4: GeV , 

for a proton or a deuteron target, and 

0.46 < Zh< 0.59 
0.68 <y< 0.86 
5.4 < VF^ < 9.3 GeV^ 

for a neutron target. The SSAs result to be almost identical to those obtained at 12 GeV; therefore, we do not show 
them explicitely. 

A word of caution is necessary when discussing JLab SIDIS observables at 6 GeV; in the high-x, relatively low-Q^ 
JLab kinematical regime, target and identified hadron mass effects (in particular for kaons), large x resummations 
and higher- twist effects might be significant. This implies potential theoretical problems in the analysis of such data. 
Notice, moreover, that a much better statistics should be expected at 12 GeV. 

V. CONCLUSIONS 

We have performed a comprehensive analysis of SIDIS data on Sivers azimuthal dependences, taking advantage of 
new and more precise experimental results. Particularly challenging are the HERMES data on kaon asymmetries, 
showing an unexpectedly large value of A^'^'^'^ "^^^ for . Our results confirm and improve previous extractions of 
the u and d Sivers distributions and offer first insights into the sea contribution to the Sivers effect. 

It turns out that the data demand a non vanishing, and large, Sivers distribution for s quarks, which, coupled to 
a new set of fragmentation functions enhancing the role of strange quarks, appears to be the only way, at present, 
to explain the data. The other sea quark (u, d, s) contributions are, at this stage, less well determined, although 
they also seem to be non vanishing. 

Taken at face value, the extracted Sivers distributions indicate a saturation of the Burkardt sum rule mainly due to 
u and d quarks alone, which carry almost opposite transverse momentum. The sea quark contribution is altogether 
rather small. This seems to rule out a contribution from gluons and, somehow, points towards a picture of the proton 
structure with the parton orbital motion restricted to valence quarks. 

Finally, we have used our extracted Sivers distributions to compute estimates for ongoing and planned new exper- 
iments at COMPASS and JLab. 
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0.02 < Ft < 1 GeV/c 
0.4 < y < 0.85 
W^> 4. GeV^ 



(27) 



0.13 <x^ < 0.40 

1.3 < Q2 < 3.1 (GeV/c)2 

2.385 <Eh< 2.404 GeV , 



(28) 
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APPENDIX A: ANALYSIS AND STATISTICAL UNCERTAINTY BANDS 



For completeness, and because it is often a matter of debate, we briefly discuss the techniques used for evahiating 
the statistical uncertainties of our estimates. A standard (see Refs. [131, (Hj) analysis is applied in order to estimate 
the values of M unknown parameters a = {ai, am}- The total is calculated by 



N 



where wc have a set of TV experimental measurements yi at known points a;,. Each measurement is supposed to be 
Gaussian distributed with variance of. The theoretical estimate F{xi]a) of the measurement j/; depends non-linearly 
on the M unknown parameters a^. 

Minimizing the total yields a set of parameters ao and a value x^j„ = x^(3-o)- 

In our particular case we have N — 173 and, for the "broken sea" ansatz, M = 11, thus ^.f. — 162; the minimum 
found by MINUIT is xfmn 162.65. 

At this stage, we would like to estimate the possible errors on the extracted parameters and the statistical uncer- 
tainties on the corresponding Sivers functions and on our estimates for the asymmetries which will be measured in 
future experiments. 

Let us call a.j some sets of parameters that could have been obtained if a slightly different set of data {xi,yi}j 
were measured (within experimental uncertainties one can generate, using Monte Carlo techniques, new data sets and 
extract the corresponding parameter sets a^). Then, for each value of j, the quantity 

Ax2^x'(a,)-x'(ao) (A2) 

is distributed according to a chi-squarc distribution, with AI degrees of freedom (see sects. 32.3.2.3 of Ref. [sj], 15.6 
of Ref. (ssj). In order to take into account the correlations between all the parameters, we would like to perform a 
joint estimation of M parameters. The corresponding coverage probability can be calculated according to the formula 

The meaning of the coverage probability can be explained as follows: suppose that we generate sets of parameters a.j, 
j = 1, ..,/, which satisfy the condition 

x'(a,)<xL„ + Ax'; (A4) 

then, if the number of sets is large enough, wc cover a hyper-volume in the M dimensional space which is called 
confidence region. The meaning of the confidence region is that with a probability P we will find the true set of 
parameters atme inside this hyper-volume (see sect. 15.6 of [1^). The minimal probability which is worth quoting is 
68.3% and is historically connected to one sigma deviation of the normal distribution, then 95.45% corresponding to 
2 sigma, etc. 

Notice that if we wanted to estimate the error of one single parameter, say oi, having fixed all the other parameter 
values, then the previous considerations would lead us back to the well known case in which M = 1 , and Ax^ = 1 for 
a required 68.3% confidence level. 

We determine the confidence hyper- volume corresponding to P = 0.9545 coverage probability for the joint estimation 
of M parameters. From Eq. (|A3[) . we obtain 

Ax^,i = 19.988. (A5) 
We then generate 200 sets of parameters a^, j = 1, 200 which satisfy the condition 

x'(aj)<xL„ + Ax/»t, (A6) 

and cover our chosen confidence region. 

Now, in all the plots the central line corresponds to the extracted set of parameters ao obtained from the xinin 
value. In order to estimate the statistical uncertainty on this result we calculate the same quantity (either single 
spin asymmetry or Sivers function) corresponding to the sets a^, j = 1, 200: at each given point x (or z or Pt, as 
appropriate) the maximal and minimal values among all of these give us the upper and lower uncertainty boundaries. 
The resulting error band corresponds to the projection of the 95.45% confidence region onto a given observable. The 
meaning of this band is straightforward: the probability to find the true result inside the shaded corridor is 95.45%. 
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The corresponding shape of the and the scatter plot of the 200 generated sets are shown in Fig. [121 where one 
can clearly see that the does not have a parabolic shape in the vicinity of ap; thus we need to take into account 
higher order corrections to its Taylor expansion. 

The explained method is general and is suitable to determine statistical uncertainties in a general case, when the 
fitting function does not depend linearly on parameters. Instead the "standard" method, exploiting the error matrix 
issued by the MINUIT package, is applicable only in the cases when the fitting function depends linearly on the 
parameters, and the parameter errors are very small. 

The statistical uncertainties in the values of {k'[) and their combinations, Eqs. have been determined 

consistently with the above procedure; these quantities have been computed for each of the 200 sets of parameters 
and the upper and lower limits shown in Eqs. simply correspond to the highest and lowest values found. 
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FIG. 9: Estimates for the single spin asymmetry Ay^j!**'* '^'^^ for pion and kaon production, which will be measured at JLab 
operating on a polarized hydrogen target with a beam energy of 12 GeV. 
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FIG. 10: Estimates for the single spin asymmetry A^^'*'' "^^^ for pion and kaon production, which will be measured at JLab 
operating on a polarized He'' (neutron) target, with a beam energy of 12 GeV. 
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target, with a beam energy of 12 GeV. 
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FIG. 12: The shape of x'^ as a function of the fit parameters (left) and the scatter plot of the 200 generated sets as a function 
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